/*
    Copyright 2006 Andrew Wilkinson <andrew@indiegigs.co.uk>

    This file is part of libfixed.

    libfixed is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    libfixed is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with libfixed; if not, write to the Free Software
    Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
*/

#include <stdio.h>

#include "fixed.h"

Fixed32 Fixed_mul32(Fixed32 arg1, Fixed32 arg2) {
    int64_t arg1m = Fixed_mag32(arg1);
    int64_t arg1f = Fixed_frac32(arg1);
    int64_t arg2m = Fixed_mag32(arg2);
    int64_t arg2f = Fixed_frac32(arg2);

    return (Fixed32) { (((arg1m * arg2m) << 32) +
           (arg1m * arg2f) + (arg1f * arg2m) +
           (((arg1f * arg2f) >> 32) & (((int64_t)1<<32)-1))) * (Fixed_sign32(arg1) * Fixed_sign32(arg2)) };
}
